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'~~> Abstract. In this note, we give constructive upper and lower bovrnds for the minimal 

iy~\ speed of propagation of traveling waves for non-local delayed reaction-diffusion 

^__( equation. 

^ 1 Introduction and the main results 

< 

fH In this note, we estimate the minimal speed of propagation of positive traveling 

'*~L wave solutions for non-local delayed reaction-diffusion equation 

I I Ut{t^ x) = Uxx{t, x) — u{t, x) + I K{x — s)g{u{t — h, s))ds, u>0, a; G M, (1.1) 



^ 



K^ which is widely used in applications, e.g. see |l2j |3j |6j |7l [91 and references 

,_^ wherein. It is assumed that the birth function g is of the monostable type, p := 

Tj- ff'(O) > 1 arid h>0. The non-negative kernel K is such that K{s) = K{—s) for 

TJ" s e M, /ju K{s)ds = 1 and /^ K{s) exp{Xs)ds is finite for all A G M. Consider 

CN 

^^ -tp^z^e) = ez'^ — z — I +pexp{—zh) / K{s)exp{—\/ezs)ds, (1.2) 

00 . . 

f^ which determines the eigenvalues of Eq. < |1.1[ | at the trivial steady state. From 

L| 1,4,2]/ we know that there is eo = eo(^) > such that ^(z, £o) = has a unique 

• ^H multiple positive root zq = za{h). Furthermore, if g{s) < g'{0)s for s > 0, then 

/\^ the minimal speed c* is equal to c* = l/^/eo ■ Note that zq and Eq are the unique 

solutions of the system 



V'(z,e) = 0, V^(z,£) = 0. (1.3) 

It is known [9] that for various systems modeled by equation | |1.1| , the minimal 
wave speed c* coincides with the spreading speed. Therefore, it is important to 
study the effects caused by the delay and other parameters (depending on spe- 
cific models) on c*, cf . f3','5"6"8 , 9 1. Another aspect of the problem concerns eas- 
ily calculable upper and lower bounds for c* . In particular, in the recent work 
IITOl , Wu et al. give several nice estimations for c* when Ka(s) = r^ — e~'' /^" 
and a < h. However, the approach of [lOJ depends heavily on the condi- 
tion a < h and on the special form of K w^hich is the fundamental solution 



of the heat equation. In the present work, we use a completely different idea 
to estimate the minimal speed for general kernels and without imposing any 
restriction on h. 

Let us state our main result. Set 



/c2 = In p / K{s) exp{ — \/lnps)ds 

Vlnp V Jr 

It is clear that fc2 > and bellow we will show that fci is positive. 

Theorem 1.1. Assume that K{s) > fs such that K{s) = K{—s) for s e M, 
J^K{s)ds = 1 and J^K{s)exp{Xs)ds is finite for all A e M. Then c* = c^,{h) = 
l/y/eo{h) is a C°°-smooth decreasing function of variable h e M+. Moreover, 



/ P - 1 2VInp f fci fca , , , , 

1. max < 2. / — — —TT -, — > < c* < mm < -, -— >, h e 0, IL 

I y p{2h+h^) + l l + h \ [l + h h ' ^ ' 



Furthermore, — < c^{h) < —, h > 1, for some positive Ci < C2. 



Observe that Theorem 1.1 implies that c^,{h) = 0(/i ^), /i -^ +00, in this 
way we improve the estimation c*{h) = 0{h'^^^), h — > +00, proved in ISllTOl. 

Proof. It follows from |T. "7] that the functions zq — zo{h) and eq = ^f){h) are 
well defined for all /i > 0. Set F{h, z, e) — (ipiz, e),'ijjz{z, e)). It is easy to see 

F e C°^iR+ xMx (0,c»),M2)^ F{h,zo,eo) = 0, and 



dF{h,zo,SQ) 



3(^0, eo) 



= ■il;zz{zo,eo)tpe{za,eo) 

= {2eo+p K{s)exp{-zo{h+ ^/eos)){h+ y/eos)'^ds) 

Jr 

-{1 + hp / i^(s) exp(-zo(/i + y/e^s))ds) > 0. 

Jr 



X ^ 

2e, 



Applying the Implicit Function Theorem we find that zq, eq £ C°°(0, +00). 

On the other hand, after introducing a new variable w = %/ez we find that 
system ||1.3| takes the following form: 



■u? I exp ( — = ] —p K(s)exp(—ws)ds. (1.4) 

/ Vve/ Jr 



w T \ ( w 



h ^ f ^ h\ I + h\ f wh 



.,2 



w^ + [2 w - 



exp ( —;= ] =p sK(s) exp(—ws)ds. (1.5) 



Let Giw) ={l+^- yj^), H{w) ={l+^-w^) exp (^) and R{w) = 
pJg^K{s)exp{—'ws)ds. Set also wo — wo{h) = ^Jea{h)zf){h). First, note that 
G{wo) = exp (^) R{wo) > and G{w) > 1 when < w < 1/Veo- As can be 
checked directly, H has a unique positive local extremum (maximum) at some 
w. Since K{s) = K{—s), s e M, it is easy to see that R increases on M.^. 



Differentiating equation | |1.4| with respect to h and using 1 1.5 1 we get the 
following differential equation 

2eoih)G{wo{h)) ^ ^ 
^"('^^ ^ 1 + hGiwoih)) > '• ^'-'^ 

The remainder of the proof will be divided in several steps. 

Step l.lihe [0, 1], then H'{1/^) = (^) e^l^° < 0. Hence, w < l/^e^. 

In addition, if w e {0,w) then H'{w) > 0. As R'{w) > for w > 0, we have 
W(i Kw < I/^/eq. Thus, we get G{wq) > 1. In this way, eo(^) > 2eo(^)/(l + h) 
for h e [0, 1] that yields (1 + hfeo{Q) < eo{h) < (1 + hfeoil)/i (equivalently, 
2c,(l)/(l + h) < c^{h) < c40)/(l + h), for h e [0, 1]). Next, taking /i = in 
equations | |1.4| and | |1.5| we obtain that 

= 2w;o(0) -p sK{s) exp {-wo{0)s)ds, (1.7) 



v^m 



l + wl{0) = p / A'(s)(l + u;o(0)s)exp(-u;o(0)s)ds 

/ Ls'^K{s)ds 2,„s fn ■s'*ivr(s)(is 4,^^ 

As a consequence of the latter formula, we get 



^o(0)< ' ^ ^ 



l + P^J^s^K{s)ds 



Then |[T7| implies that c*(0) < h so that c*(/i) < ki/{l + h) for h < I. Note 
that /ci > since R is increasing for it; > 0. Finally, since c* (/i) is decreasing, 
we have that c^{h) < fci/2 for h> 1. 

Step II. If /i > 1, then ui > l/^eo. As consequence, G('u)) < 1 — G{l/y/eo) 
so that G'(fi;) > G{w) for all w e [0,w] (see Figure [ij. Additionally, G(w) = 
{2wy/eo — 1)^ > ^, therefore we conclude that G{wo) > l/h. Hence, we have 
s'oih) > eo{h)/h, so that e{h) > e{l)h (equivalently, c*{h) < c*(l)/v^) for h > 



1. Now, if ft, = 1 we have w = l/\/eo(l)- Thus, taking h = 1 and w = w in ( 1.4 



wegetexp(l/eo(l)) = -R(w) > -R(O) =p that yields ^75^ < l/^eo(l) = c*(l). 
On the other hand, for all < it; < \/ y/e^, we have 

exp ( -^= ) < ( IH ^ - w^ j exp ( -^= ] <P j K{s) exp(-ws)(is. (1.8) 

1 and w = 



In particular, taking h = 1 and tu = ^/Inp in 1 1.8 1 we conclude that c* (1) < fc2 
so that c^,{h) < k2/VhioT h > I. Additionally, using c*(/i) > 2c*(l)/(l + h) 
obtained in step I, we also concluded that c^{h) > 2^/]iip/{l + h), for h e [0, 1]. 




Figure 1: G, H and R for h>l. 

Step III. For /i > 0, it is evident that egC^) < 2eo{h)/h. Integrating the latter 
inequality on [h, 1] we obtain e{h) > e(l)/i^ (equivalently, c^{h) < c*(l)//i), for 
< h < lsothatc*(/i) < k2/h,iorh e (0,1]. Analogous, by integrating eo(/i) < 
2eo{h)/h on [1, /i] we have eo(/i) < eo{l)h'^ (equivalently, c^,{h) > c*(l)//i), for 
h> I. Thus, we obtain c* (h) > \/\np/h, h>\. 

On the other hand, for all /i > 0, we have G(wo) < 1 + l/(4eo)- As 
consequence, eo(/i) < {Aeo{h) + 1)/(2(1 + /i)) for all /i > so that eo(/i) < 
((4eo(0) + 1)(1 + hf - l)/4. Taking h = in |L4) , we get 1 + l/(4eo(0)) > 
G(wo(0)) = i?(wo(0)) > _p so that c*(0) > 2Vp^^l7Tn consequence. 



c*(/i) > 2, 



p-l 



p{2h + h'^) + l' 



h>0. 



(1.9) 



Step IV. Setting w = r,r £ (0,1), in the second inequality of 1 1.8 i we obtain 



(1 - r ) exp 

from which we get that 

1 1 , 

v/i^ hr VI 



rh 



V^oW) 



<p K{s)exp{—rs)ds, 

Js. 



K{s)exjp{—rs)ds], h>0 



(1.10) 



I I I c c 

Considering ( 1.9 1 and 1 1.10 1 we get — < c^,{h) < — ior h> 1. This completes 

the proof. D 



2 An example 

Consider the heat kernel if Q,(s) = (47ra)~^/^ exp (— s^/(4q;)). Then Theorem 
1 1.1 1 applies with 



ki^2^p-l 



I + ap exp 



/a(p-i)y 



(1 + a)\/lnp. 



a/1 + ap 
In fact, in this case we can plot graphs of c* against h using standard numerical 
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Figure 2: The minimal speed and its bounds {p — 2 and a = 1). 

methods to solve some appropriately chosen initial value problem eo(^o) — Po 
for differential equation ( 1.6| . For example, if we take ho — a then po coincides 
with positive solution of the equation 1 + -p = pexp{—f-). Next, we can 



explicitly find G{wq) in 1 1.6 1 by using Cardano's formulas to solve the cubic 
equation {wq — wo/^/eo — l){2y/eoawo — h) + 1 — 2^eowo ~ 0. It is easy to see 
that this equation has three real roots for all /i > and a > 0, and that wq is the 
leftmost positive root. 

Figure [2] shows the minimal speed c* and its estimations when p = 2 and 
a = I. Remark that we do not need the restriction a < h required in IITOll . 

Finally, note that letting a -^ 0+ in 1 1.1 1 and 1 1.2 i we recover the character- 



istic equation for the delayed reaction-diffusion equation 

ut{t, x) = Uxx{t, x) - u{t, x) + g{u{t - h, x)), 



which was studied by various authors (e.g. see fV.'E] and references therein). 
In this case, our results complete and partially improve the estimations of |j8l. 
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